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In many sliding systems consisting of solid object on a solid substrate under dry condition, the friction force 
does not depend on the apparent contact area and is proportional to the loading force. This behaviour is 
called Amontons' law and indicates that the friction coefficient, or the ratio of the friction force to the 
loading force, is constant. Here, however, using numerical and analytical methods, we show that Amontons' 
law breaks down systematically under certain conditions for an elastic object experiencing a friction force 
that locally obeys Amontons' law. The macroscopic static friction coefficient, which corresponds to the onset 
of bulk sliding of the object, decreases as pressure or system length increases. This decrease results from 
precursor slips before the onset of bulk sliding, and is consistent with the results of certain previous 
experiments. The mechanisms for these behaviours are clarified. These results will provide new insight into 
controlling friction. 

When we apply a shear force to a solid object on a solid substrate to start a sliding motion, the shear force 
must be greater than the maximum static friction force. When the object is sliding, the kinetic friction 
force applies. Friction plays an important role in various phenomena ranging from those at the 
nanometre scale to earthquakes; the phenomenon of friction has been investigated since ancient times 1 " 5 . 
From the engineering point of view, friction is required to be small in certain cases and large in other cases. 
The control of friction is one of the key factors towards achieving improvements in green technology and 
nanotechnology. In the 15th century da Vinci discovered that the friction force is proportional to the applied 
loading force and is independent of the apparent contact area between two solid surfaces 1 " 5 . This behaviour of 
friction was rediscovered by Amontons approximately 200 years later; this law is now called Amontons' law of 
friction and holds for various systems at first approximation 1 " 5 . The ratio of the friction force to the loading force 
is called the friction coefficient, and according to Amontons' law, this coefficient does not depend on the loading 
force or the apparent contact area. The mechanism of friction was explained in the mid-twentieth century by 
Bowden and Tabor 1 . For actual solid surfaces in contact with each other, because of surface roughness, only a tiny 
fraction of the surfaces form junctions, the so-called real contact points. Amontons' law is explained as resulting 
from the increase in the total area of real contact points, that is, the real contact area, in proportion to the loading 
force and the constant binding energy per unit real contact area 1 " 7 . However, the mechanism and the validity of 
Amontons' law are still discussed actively 2 " 13 . 

Another interesting problem concerning sliding friction is the question of how a macroscopic object begins to 
slide. Usually, it is considered that a shear force smaller than the maximum static friction force does not induce 
any slip motions. However, recent measurements of the instantaneous local real contact area density of poly- 
(methyl methacrylate) (PMMA) show that precursors appear as local slips at the interface under shear forces well 
below the maximum static friction force, and this force corresponds to the onset of bulk sliding 14 " 17 . If the shear 
force is applied slowly from the trailing edge of the sample, a discrete propagation sequence for the local slip 
appears 15 . Each front of the slip starts from the trailing edge and stops after propagating a certain length greater 
than that of the previous front. When the slip front reaches the leading edge, bulk sliding occurs. Similar 
behaviour is observed in numerical studies based on ID 17 " 20 and 2D 21 spring-block models. The variation in 
the front velocity and the front velocity's dependence on the stress distribution have been observed in experi- 
ments 1416 and have been examined using numerical approaches 18 ' 21 ' 22 . A slow precursor with a finite front velocity 
has been investigated with a ID continuum model 23,24 . The precursor dynamics are expected to relate to the 
inhomogeneity of the system and are expected to be important to the understanding of friction. In fact, ref. 9 
reported that the macroscopic static friction coefficient /n M , which corresponds to the maximum static friction 
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force, varies with the experimental loading configuration, and this 
variation is related to the precursor dynamics. Similar behaviour for 
the maximum static friction force has also been investigated with a 
ID spring-block model 25 . Despite these studies, however, the mech- 
anism of the precursor and its relation to the maximum static friction 
force have not been clarified theoretically. 

In this study, using both numerical and analytical methods, we 
investigate the precursor dynamics and the relation of these 
dynamics to the macroscopic static friction coefficient fi M of an 
elastic object in contact with a rigid substrate, and we show that 
/n M decreases as pressure or system length increases. This behaviour 
indicates the systematic breakdown of Amontons' law. In the present 
system, the elastic object is subject to viscous damping and a friction 
force that obeys Amontons' law locally; however, the law breaks 
down for the overall system. The breakdown of the law results from 
a quasi-static precursor appearing in the propagation of local slips 
before the onset of bulk sliding and from the transition of the quasi- 
static precursor to a dynamic precursor at a certain critical precursor 
length. For sufficiently large or small pressures or system lengths, fi M 
becomes nearly constant and Amontons' law approximately holds. 
The relation of /n M to the critical precursor length and the mechanism 
of the breakdown are clarified analytically with a ID effective model. 
These behaviours arise from the distribution of pressure on the bot- 
tom of the object, which results from the torque induced by the shear 
force, and from the competition between the viscous damping and 
the velocity- weakening friction. Consequently, the behaviours do not 
depend on the details of the system. It is known that Amontons' law 
does not apply under some conditions, such as in the presence of 
strong adhesion, in the absence of multiple real contact points, with 
nonlinear dependence of the real contact area on the applied load, 
with surface detachment, with a change in surface conditions, and 
for gels 2 " 4 ' 61011 . However, the mechanism of the breakdown of 
Amontons' law discussed in this paper is different from those known 
previously. In addition, a theoretical prediction of the behaviour of 
the friction coefficient is given for the first time. These results will 
provide new approaches to control friction. 



Results 

Finite element method (FEM) calculation. We first analyse the 
friction behaviour for an elastic block with length L, width W, and 
height H along the x-, y-, and z-axes, respectively, on a rigid substrate 
(see Fig. 1(a)) using a 2D finite element method (FEM). We assume 
that the block has viscous damping proportional to the strain rate 
and has friction stress proportional to the pressure at the interface; 
that is, the local friction force obeys Amontons' law. Note that the 
resolution of the models employed in this work is larger than the 
mean distance between the real contact points between the two 
surfaces and smaller than the length scale of the variation of the 
stresses discussed below. The validity of the law on this length 
scale for PMMA is supported by the proportionality of the real 
contact area to the loading pressure 7 , with the assumption of a 
constant binding energy per unit real contact area, and the law's 
validity is consistent with the friction experiments in which the 
size of the apparent contact area is on the order of 1 mm 2 in a 
multiple real contact point configuration 6 . The local friction 
coefficient fi (u x ) is given by the local static friction coefficient jd s 
when the local slip velocity u x is equal to zero, and it decreases 
linearly from [i s to the local kinetic friction coefficient fi K with 
increasing u x until the local friction coefficient finally equals ^ K at 
velocities greater than characteristic velocity v c . The plane stress 
condition is assumed along the y-axis. We apply a uniform 
pressure P ext = F N /(LW) to the top surface, where F N denotes the 
loading force. The shear force F r is applied from the trailing edge at 
height h by means of a rigid rod with a width of 0.1H. The rod moves 
with a constant velocity V that is sufficiently slow. We set W = 1 and 
h = HI 2. The results shown here do not depend on the details of the 
system parameters. However, for the quantitative calculations, we 
employ two parameter sets; parameter set (A) simulates a model vis- 
coelastic material, and (B) simulates PMMA. Hereafter, we employ 
quantities normalised by the mass density, Young's modulus, and 
viscosity, and these quantities are expressed with a tilde. The 
following results are obtained for d = H /Z = 0.5. (See Methods for 
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Figure 1 | (a) An elastic object under uniform pressure P ext = F N /(LW) on a rigid substrate is pushed at height h by a rigid rod with shear force F T . 
(b) The ratio Ft /F N as a function of U = Vt. The horizontal line indicates the value of fi M . The inset shows an enlarged view of the box indicated by the red 
outlines, (c) The local slip region is shown in red in the U—x plane. The horizontal line indicates the critical length of the quasi-static precursor £ c . (d) The 
normalised instantaneous local density of the real contact area in the region of the box indicated by blue outlines in (c) . The thin line indicates the position 
of the precursor front. The results are obtained for L= 1.0, P ext = 0.003 and parameter set (A). 
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more details regarding the FEM calculation, the velocity dependence 
of the local friction coefficient, and the parameters.). 

Figure 1(b) shows the ratio F T /F N as a function of the displace- 
ment of the rod U = Vt, where t denotes time. The ratio repeatedly 
shows a nearly linear increase and a periodic large drop after the first 
drop. In this work, we focus on the periodic regime. The periodic 
behaviour corresponds to a periodic stick- slip motion. The large 
drop is caused by a large slip accompanied by bulk sliding. The value 
of F T at the peak of the F T / Fn curve is the maximum static friction 
force. The inset shows an enlarged view of the box indicated by the 
red outline. A sequence of small drops in F T j F N is clearly observed 
before the onset of bulk sliding at U~ 0.0047. As discussed below, 
each small drop is accompanied by a rapid local slip. This behaviour 
has previously been observed in certain experiments 15 and numerical 
studies 17 " 21 . Figure 1(c) shows the local slip region, where the slip 
velocity at the interface has a finite value, in the U—x plane. Here, x 
denotes the position at the bottom of the object. At values well below 
the maximum static friction force, the shear force Fr induces a slow 
precursor slip from the trailing edge, and this precursor slip is called 
the quasi-static precursor. The velocities of the local slip and the front 
of this precursor are proportional to the driving velocity, and they are 
vanishingly small. The precursor length £ increases with increasing 
U. When £ reaches a critical length £ c , the quasi- static precursor 
becomes unstable and transforms to a leading rapid precursor with 
large slip and front velocities. The leading rapid precursor nucleates 
near the trailing edge. The front of the leading rapid precursor begins 
with a velocity close to that of sound and in a manner similar to the 
supershear rupture observed in experiments 916 . Subsequently, the 
front velocity reduces and approaches the Rayleigh wave velocity. 
When the front enters the leading edge, bulk sliding occurs, and F T 
shows a large drop in value. 

Figure 1(d) shows the normalised instantaneous local density of 
the real contact area in the region of the box indicated in (c). The local 
density is calculated from the local slip velocity and decreases with 
increasing velocity. (See Methods for this calculation.). Before the 
onset of bulk sliding at U~ 0.0047, a discrete sequence of rapid 
precursors, similar to that observed in experiments 15 , appears at 
U~ 0.0045, and this set of precursors causes the sequence of small 
drops in Ft / Fn observed in the inset of (b). This type of precursor is 
called the bounded rapid precursor. Each precursor nucleates in the 
region x < £ with a front velocity that is close to the velocity of sound 
and is independent of the driving velocity and nucleates in a manner 
similar to the leading rapid precursor. Subsequently, the front decele- 
rates and stops after propagating a certain length. The propagation 
length of the precursor front increases with U. This increase is also 
observed in experiments 15 . The local slip velocity of this precursor 
is considerably larger than that of the quasi-static precursor shown 
in (c); however, this slip velocity decreases with decreasing driving 
velocity in contrast to the behaviour of the slip velocity for the leading 
rapid precursor. The quasistatic precursor observed in (c) almost 
disappears for the real contact area density shown in (d) because 
of its vanishingly small slip velocity. The disappearance of the 
quasi-static precursor is also consistent with the results of certain 
experiments 15 . 

We define the macroscopic static friction coefficient [i M to be the 
peak value of F T /F N in Fig. 1(b). The pressure (Pext) dependence of 
jd M is shown for various values of the system length L in Figs. 2(a) and 
2(b) for parameter sets (A) and (B), respectively. Clear decreases in 
fi M are observed with increasing P e xt or L values. It can be shown that 
jd M also depends on the apparent contact area. These behaviours 
indicate the breakdown of Amontons' law and the extensive property 
of the friction force. This pressure dependence is consistent with the 
results of experiments conducted with PMMA 8,9 . The magnitudes of 
the normalised parameters, P ex t and L, depend on the mass density, 
Young's modulus, and the viscosity coefficient. Therefore, in contrast 
to the general belief that the friction coefficient depends mainly on 



the surface properties, fi M also depends on these bulk material para- 
meters, as noted in ref. 25. The variation in fi M results from the 
formation of the precursors before the onset of bulk sliding. 
Figure 3(a) shows that fi M is scaled by the critical length of the 
quasi-static precursor normalised by the system length, £ C /L, for 
each parameter set. The dependence of £ C /L on P ex t> shown in the 
inset, with this scaling of [i M on £ C /L indicates the P e xt dependence of 
Mm- 

Figure 3(b) shows the pressure at the interface p(x) for the mag- 
nitudes of U indicated by arrows in Figs. 1 (b, c). The applied pressure 
at the top surface is uniform; however, p(x) is not uniform and 
increases with x. Similar pressure distributions are observed in 
experiments 9,16 and result from the torque induced by the shear force 
F T 19 . The low pressure causes the local maximum static friction force 
to be small. As a result, the quasi-static precursor starts from the 
trailing edge at x = 0. Figure 3(c) shows the ratio of the shear stress at 
the interface t(x) to the pressure, t(x) /p(x). Immediately after bulk 
sliding stops, this ratio approximately equals the value of ^ K (dashed 
line) for the entire interface because of the large local slip velocity 
accompanied by the bulk sliding and the finite relaxation time of the 
local stress. The magnitude of t(x)is equal to the local friction stress 
MkP(^) at me instant of the vanishing acceleration of the rapid local 
slip, and this magnitude has almost no change during the decelera- 
tion and after stoppage because of the finite relaxation time. 
However, in the region of length £ where the quasi -static precursor 
front has passed, z(x)/p(x) equals /n s (straight line) because t(x) is 
given by the local friction stress for vanishing velocity ^siK*)- When 
£ reaches the critical length £ c , the quasi-static precursor transforms 
to the leading rapid precursor. Subsequently, the front of this rapid 
precursor enters the leading edge of the system quickly and bulk 
sliding occurs. The stress distribution has almost no change during 
the propagation of this front because of the short duration of the 
propagation. As a result, fi M is determined by £ c , and this relation is 
shown analytically below. Note that the stress relaxes slightly after 
the appearance of each of the bounded rapid precursors, and the 
stress recovers its original value quickly through the following 
quasi-static precursor, as shown in Fig. 1(d). Figure 3(c) also shows 
that T(x)/p(x) can be larger than the macroscopic static friction 
coefficient /n M without precipitating any local slip, as was observed 
in experiments 916 . 

Analysis based on a ID effective model. To analyse the abovemen- 
tioned numerical results, we employ a ID effective model in which 
the degrees of freedom of the elastic object along the z-axis are 
neglected. The equation of motion of the model is expressed by 




0 0.02 0.04 0 1.5x 10 



± ext 1 ext 

Figure 2 | Macroscopic static friction coefficient /i M as a function of P ex t 
for parameter sets (A) (a) and (B) (b). The lines with symbols indicate the 
results of the FEM calculation. Thin lines indicate analytical results based 
on the ID effective model, where we set a = 0.2. Lines of the same colour 
correspond to the same value of L. 
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Figure 3 | (a) Macroscopic static friction coefficient fi M as a function of £ C /L for various values of P ext and L. The open symbols indicate results for the 
parameter set (A) with L = 0.25 ( □ ), 1.0 ( O ), and 4.0 ( A) and the filled symbols indicate those for (B) with L = 0.02 ( ■ ), 0.04 (#), and 0.08 (A). The 
lines show the theoretical result as obtained with equation (2). The inset shows the P ex t dependence of £ C /L for parameter set (A), (b) The normalised 
pressure /P e xt and (c) the ratio of the shear stress to the pressure x(x) jp{x) at the interface for the magnitudes of U indicated by arrows in Figs. 1(b) 
and 1 (c). The arrows in (c) indicate the positions of the precursor front. The three horizontal lines ji s , /x M , and ji K in order from the top of the panel. The 
parameters are the same as those for Figs, l(b-d). 



d 2 u(x~t) _ dd xx (x~t) a xz (x~t) - fi(u)p(x) 



dt 2 



dx 



oiH 



Here, u(x~t) denotes the displacement along the x-axis at the inter- 
face, fi(u)p(x) is the friction stress, and ocH denotes the characteristic 
length of the variation in the (xz) component of the stress. We 
employ a as a fitting parameter. The normal and the shear stresses 
are respectively given by a xx =Eidu/ dx + ij t d 2 u/ dxdt, and a xz = 

E 2 (p — u) + (V — du/dt) 1 2h y where rj t denotes the effective vis- 
cosity and Ei and E 2 represent effective elastic constants (see 
Methods for these parameters and the boundary condition.). We 
assume that the pressure at the interface is given by 
p(x) = 2P ext x/L, which simulates the FEM result shown in 

Fig. 3(b). Hereafter, we set the origin of U to be the position of the 
pushing rod just after bulk sliding stops. The adiabatic solution of 
equation (1) with a precursor of length £ is obtained analytically, and 
the solution gives £~E 2 LU ^ ^(fi s — ^ K )/zP ex t} for £/L<^l (see 

Supplementary Information). Hence, £ increases adiabatically with 
an adiabatic increase in U. This adiabatic increase in length indicates 
the precursor is quasi-static. A similar relation between the precursor 
length and the shear force is obtained with a ID spring-block model; 
however, in this case the precursor is not adiabatic 20 . When £ reaches 
a critical length £ c , the quasi-static precursor becomes unstable and 
transforms to the leading rapid precursor. The leading rapid pre- 
cursor leads to bulk sliding. Substituting the adiabatic equation and 
the relation between U and £ into the expression for the shear force 
= dxa xz (x) and setting £ = £ c ,we obtain 



0m = /*k + O*s-/*k)4/^ 



(2) 



for £ C /L<^1. As shown in Fig. 3(a), this relation agrees well with the 
FEM calculation, even for £ C /L< 1. 

The critical length £ c is obtained by a linear stability analysis of 
equation (1) in the limit of vanishing V (see Supplementary 
Information). The analysis provides the equation for the n-th eigen- 
value cb n of the time evolution operator for the fluctuation by 



1 r 



v c ocd 



L I 2\L 

2 

1 (r 



k z \L 



P ext L ( ~ ]cb n = 0, 



(3) 



where k 2 =Ei/ E 2 , r 2 = ochii, and k n = (n + 1/2)tl Here, the off-diag- 
onal terms, which are of a higher order in £, are neglected. The 
eigenvalue w n gives the instability condition of the adiabatic solution. 
The instability results from the appearance of a positive real part for 
any eigenvalue. For small values of £/L, the viscosity expressed by the 
second term in equation (3) stabilises the adiabatic motion. However, 
the velocity- weakening friction force expressed by the last term leads 
to instability in motion for large £/L. In the present case, the oscil- 
lating fluctuation corresponding to a complex w n with a positive real 
part never grows, because the backward motion of the oscillation 
relaxes the local shear stress, and this motion is inhibited by the 
friction stress. Instead, this oscillating fluctuation yields the bounded 
rapid precursors observed in Fig. 1(d). 

The appearance of the positive real part of (b 0 at £ = £ sc , where £ sc 
denotes the subcritical length, yields the first bounded rapid pre- 
cursor, and the appearance of cb n with n > 1 is considered to yield 
each of the subsequent precursors appearing at £>£ sc .^ Further 
increase in £ causes w 0 to become real and positive at 4>4o a 
positive, real w 0 results in the growing instability and causes the 
leading rapid precursor and the subsequent bulk sliding. The relation 
between £ c and P ex t is given by the instability condition, (b 0 > 0. As 
mentioned above, the instability is caused by the competition 
between the viscosity and the velocity-weakening friction force, 
which are expressed by the second and last terms in equation (3), 
respectively. Hence, £ C = L decreases as P e xt or L increases because the 
last term is enhanced by P ext L; for £ C /L<^i 1 , £ C /L oc (LP ext ) ^ 3 as is 
easily seen from equation (3) (also see Supplementary Information). 
This result is consistent with that of the FEM calculation shown in 
the inset of Fig. 3(a). By inserting £ c (P ext ,Vj, which is given by the 
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abovementioned relation, into equation (2), we obtain the pressure 
and system length dependence of /j M , and /n M decreases with increas- 
ing P ext or L. The results are shown in Figs. 2(a,b). The analytical 
results agree with the FEM calculation semiquantitatively for para- 
meter set (A) and qualitatively for parameter set (B). The deviation of 
the analytical results from the FEM calculation may arise from the 
absence of internal degrees of freedom along the z-axis in the ID 
effective model given by equation (1). 

Discussion 

In this work we observed three types of precursors prior to the 
occurrence of bulk sliding: the quasi-static precursor and the 
bounded and the leading rapid precursors. The latter two are seen 
to correspond to those observed in experiments 914 " 17 . The quasi- 
static precursor is also observed in numerical studies 23 , but not in 
experiments 914 " 17 , because it vanishes away in the local density of the 
real contact area measured in experiments, as previously discussed. 
Some numerical studies 1719 " 21 do not report this precursor. The lack 
of this precursor is because local friction coefficient discontinuously 
decreases as velocity increases in these studies. This discontinuous 
decrease in the local friction coefficient corresponds to a vanishing v c 
and a subsequently vanishing stable region for the quasi- static pre- 
cursor. In ref. 22, a model similar to the present FEM model is 
employed; however, the value of the viscosity coefficient corresponds 
to £ C /L~0. This result is also consistent with the absence of the 
quasi-static precursor. 

In ref. 9, it was reported that ^ M varies with a certain acceleration 
length of the leading rapid precursor. The acceleration length is 
difficult to define in the present work, and this difficulty may result 
from a loading configuration different from that of the experiment. 
The relation of ^ M with the acceleration length will be studied in 
future work. However, the measured values of /n M in ref. 9 show a 
tendency to decrease with an increase in the loading force and 
depend on the stress distribution in the system. These behaviours 
are consistent with those observed in the present work. 

Earthquakes are among the largest friction phenomena on the 
earth. Many large earthquakes are preceded by foreshocks. In the 
2011 Tohoku-Oki earthquake with magnitude 9, the propagation of 
slow slip preceding the main shock was observed over approximately 
1 month 26 . The behaviour is similar to the propagation of the pre- 
cursor front observed in this work. The present results will provide 
new insights into earthquake behaviour. 

In conclusion, using both numerical and analytical methods, we 
show that the macroscopic static friction coefficient ^ M of an elastic 
object decreases with increasing pressure or system length. The mag- 
nitude of /n M also depends on the apparent contact area and the bulk 
material parameters. These behaviours indicate the systematic break- 
down of Amontons' law. The elastic object is subject to viscous 
damping and a friction force that obeys Amontons' law locally; how- 
ever, the law undergoes breakdown as a whole. The behaviour of fi M 
is consistent with the results of relevant experiments 8,9 , and it arises 
from the occurrence of precursor slips before bulk sliding and from 
the transition from quasi- static to rapid motion of the precursors. 
The linear stability analysis based on a ID effective model gives the 
relation between fi M and the critical length of the transition, and it 
clarifies the mechanism of the transition. The transition is caused by 
the pressure distribution at the bottom of the object, which results 
from the torque induced by the shear force, and by the competition 
between the viscous damping and the velocity- weakening friction. If 
the ratio of the critical length of the transition to the system length is 
considerably less than or approximately equal to unity, ^ M becomes 
almost constant and Amontons' law holds approximately. These 
situations are caused by sufficiently large or small values of pressure 
or system length. The qualitative features of these results do not 
depend on the details of the model for the case in which the pressure 
at the bottom of the object increases monotonically along the driving 



direction, the local friction coefficient decreases with velocity con- 
tinuously, and the system has viscous damping. According to these 
results, an object with greater width or an object in which the part 
contacting with the interface is divided into shorter portions along 
the driving direction can have a larger maximum static friction force 
under a constant loading force. The present results will provide new 
techniques for controlling friction. 

Methods 

FEM calculation. The equation of motion of the displacement vector of the elastic 
object u(r, t) employed in the FEM calculation is expressed as pd 2 u(r,t) j dt 2 = V-a, 
where p is the mass density. The stress tensor a is composed of an elastic part <7* el) and 
a dissipative part. The dissipative component is expressed as 
ff i/f =y li£xP + r l2(cxx + €zz)daiJ- Here e aj 6 denotes a component of the strain rate 
tensor, and rj ij2 denote the viscosity coefficients 27 . The boundary condition is given by 
g zz = —P ex t and o xz = 0 at the top surface, and o xx = 0 and a zx = 0 at x = 0 and L 
except for the portion in contact with the rigid rod. 

In the FEM calculation, the object is divided into equal- sized rectangular cells. The 
number of cells is 40 X 40 or 80 X 80. The convergence of the results with respect to 
the FEM mesh size is verified. At the beginning of the FEM simulation, we apply a 
uniform pressure to the top surface causing the elastic object to relax. Subsequently, 
we start pushing the object at t = 0. The results shown in the paper are obtained for 
V = 10~ 5 . This value of V corresponds to the adiabatic limit except for the results 
shown in the inset of Fig. 1 (b) and in Fig. 1 (d). 

Velocity dependence of the local friction coefficient and the real contact area 
density. The velocity dependence of the local friction coefficient employed in this 
work is obtained from the rate- and state-dependent friction law, \i = fi K + (fi s — 
li K )(j) with 0 = (1 — 0)/t— \u x \/D, as discussed in ref. 4. Here, 0 < (j) < 1 represents 
the state variable and x and D denote the relaxation time and length, respectively. The 
above equation yields fi(u x ), as noted in the text for v c = D/t, in the limit of vanishing 
t, where (j) = l — \u x \/v c . The magnitude of (f) is proportional to the deviation of the 
local density of the real contact area from that for | u x \ > v c . Thus, we obtain the 
normalised instantaneous local density of the real contact area shown in Fig. 1(d). 

Parameters. We introduce the normalised quantit ies in which the length, mass and 
time are normalised by L 0 = *h / y/pEj = q\ / \J pE 3 and to = rj 1 /E, respectively, 
where E represents Young's modulus. In this paper, the normalised quantities are 
expressed with a tilde. We employ two parameter sets: parameter set (A) simulates a 
model viscoelastic material, which has Poisson's ratio v = 0.34, fi s = 0.38, fi K = 0.1, 
v c = 3.4 x 10 " 4 , and rj 2 = l; and (B) simulates PMMA, which has v = 0.4, ^ s = 1-2, 
jU K = 0-2, v c = 3.9 x 10~ 7 , and rj 2 = 1- The present FEM calculation shows that /u s and 
fi K correspond to the maximum and minimum values of the ratio of the local shear 
stress to the pressure, respectively. Hence we estimate ^ s and fi K for PMMA from the 
maximum value of the ratio in Fig. 4(c) and the minimum value of the ratio in 
Fig. 4(b) in ref. 9, respectively. The value of v c is estimated from the values of fi s , A*k 
and the linear fit of the steady- state velocity dependence of the rate- and state- 
dependent friction law reported in ref. 5. The value of L 0 for PMMA is estimated as 
follows. In ref. 9, ^ M = 0.6 for L = 0.2 m, d = 0.5, andP ext = 1.67 X 10 6 Pa. The FEM 
calculation yields fi M = 0.6 for P ext L = 2.3 x 10 ~ 6 . Consequently, we obtain L 0 = 
48 m. These values of the parameters for PMMA have some uncertainties. However, 
the essential feature of the results obtained here is, as mentioned before, not specific to 
these values. The local validity of Amontons' law for PMMA may require some 
discussions. Even if Amontons' law breaks down locally, the behaviour of the 
precursor dynamics and the static friction force discussed here will be still 
qualitatively valid because their mechanism is general. For metals, Amontons' law is 
expected to better hold locally and the present qualitative results will also be 
applicable. 

ID effective model. In the ID effective model, the effective viscosity is given by 
fj t = 1 + r\ 2 and the two elastic constants are given by Ej = l/{ ( 1 + v) (1 — v) } and 
E 2 = 1 / {2(1 + v) }. The boundary condition is given by du(x,t) J dx = 0 at x = 0 and L. 
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